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Abstract 

Signal of the rare radiative decay 4> — > r]ir°j 
was observed in the SND experiment at 
VEPP-2M electron-positron collider. The re- 
sult is based on the analysis of data, corre- 
sponding to a total integrated luminosity of 
4po _1 , or 8 • 10 6 ^-mesons produced. The 
measured branching ratio of cf) — > 7]ir°j decay 
is equal to (0.83 ± 0.23) • 10~ 4 . 

1 Introduction 

Several mechanisms may contribute into 
0(1020) — > 777r°7 decay rate including cf) — > 
Vw° type transitions, where V is a vector 
meson, with a subsequent V — > iff radiative 
decay, and electric dipole radiative transition 
<fi — ► ao (980)7 with a scalar ao(980) meson 
decaying into tjtt . 

Theoretical predictions of the <fi — > rjir 7 
branching ratio, based on assumption of pure 
cj> -> Vir° decay mechanism (|l), §, f§), are 
quite low, of the order of 5 • 10 -6 . Situa- 
tion with theoretical description of the (f> — * 
ao (980)7 decay is much more complicated, 
because the results strongly depend on the 
quark structure of the ao(980) meson. At 
present, its quark structure is not well es- 
tablished and several models exist including 
modifications of qq scheme Q , 4-quark model 
||, and KK molecular model [||. Theoreti- 
cal estimations of cf) — > ao (980)7 branching 
ratios vary from about 10~ 5 in two-quark 
and KK molecular models up to 10~ 4 in 
the 4-quark model ( , @] , fl ) , making the 
<f> — » ao (980)7 decay a unique probe of the 
structure of ao(980) meson. 

The first search for the 4> — * r)ir°j decay 
was conducted in the ND experiment || at 
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VEPP-2M e+e" collider in 1987, where an 
upper limit of 2.5 • 10~ 3 was established. In 
1995 the new SND detector having better 
hcrmeticity, granularity, energy and spatial 
resolution started operation at VEPP-2M. 
First indications of </> — > r)TT 0r y decay were 
seen by SND in 1997 @. The result was 
based on analysis of half of the 0-meson data 
recorded in 1996-1997. In this work analy- 
sis of full data sample was carried out. The 
4> — > ry7r°7 decay was observed for the first 
time and its branching ratio was measured. 

2 Detector and experi- 
ment 

Experiment (To), |ll| was carried out in 1996 
at VEPP-2M e+e" collider with SND detec- 
tor. The SND detector jl2| is a universal 
nonmagnetic detector. Its main part is a 3- 
layer electromagnetic calorimeter consisted of 
1630 Nal(Tl) crystals with a total thickness 
of 35 cm or 14 radiation lengths. The en- 
ergy resolution for photons can be described 
as 4.2%/ i y /E{GeV) fl§, the angular reso- 
lution of the calorimeter is close to 1.5 de- 
grees, and the solid angle coverage is 90% of 
47r steradian. 

The data, used for the study of (f> ~ + r/TT ^ 
decay was recorded in 1996-1997 in the en- 
ergy range from 980 to 1040 MeV with the 
most of the data collected in the close vicinity 
of the 0-peak. Six successive scans of the en- 
ergy range were performed with a total inte- 
grated luminosity of A.Opb^ 1 and total num- 
ber of 0- mesons produced of 8.2 • 10 6 . 



3 Data analysis 

The main background sources for the process 
under study 

e + e~ — * (f> — » i]tt°j — > 57 (1) 
are the following 0-meson decays: 

e + e~ — > — ► 777 — * 37r°7 (2) 

e + e~ — > — > KsKl — ► neutrals (3) 
e + e~ -> -> 7r°7r°7 (4) 
and nonresonant process 

e+e- W7T° 7r 7r° 7 . (5) 

The expected number of events of the pro- 
cess (1) at a branching ratio of <fi — > 7771-° 7 
equal to 10 -4 is about 300, while the number 
of background events (2) is 3 • 10 4 . Although 
the process (2) does not produce 57 events di- 
rectly, it can fake the topology of the process 
(1) due to either merging of close photons or 
loss of soft photons through openings in the 
calorimeter. The process (3) produces 8 • 10 5 
KsKl events with Ks — * 7r°7r° decays. The 
Kl-s in </>-meson decays are slow and have 
a decay length of 3 m, while the nuclear in- 
teraction length in Nal(Tl) is about 30 cm. 
Characteristic feature of spurious 57 events 
produced by K S K L decays due to either nu- 
clear interaction of Kl-s or their decays in 
flight is an energy-momentum imbalance and 
poor quality of at least one photon in the 
event. 

Primary event selection was based on sim- 
ple criteria: the number of reconstructed 
photons is equal to five, there are no tracks 
in the central drift chambers, the total en- 
ergy deposition ranges from 0.8 to 1.1 cen- 
ter of mass energy 2Eq, and the total trans- 
verse momentum is less than 0.15Eq/c. Such 
criteria greatly reduce background from the 
processes (2) and (3), not affecting the de- 
cays (1), (4), and (5). Also rejected were 
the events containing close photon pairs with 
spatial angle between photons less than 27 
degrees, where energies could be poorly re- 
constructed. 

Next step in the event selection was based 
on quality of photons and kinematic fit. Pa- 
rameter ( was constructed describing the 
likelihood of a hypothesis, that given trans- 
verse energy deposition profile of a cluster 
of hit crystals in the calorimeter can be at- 
tributed to a single photon. The require- 
ments were imposed that £ < 0, for all pho- 
tons. Then kinematic fit was performed un- 
der the assumption, that selected events are 



e + e~ — > 57 ones and corresponding param- 
eter Xei describing the likelihood of this as- 
sumption was calculated. Events with xl > 
10 were also rejected. Study of simulated 
57 events of the processes (1), (4), and (5) 
showed, that the \\ anQl C cut s reject less 
than 15% of true 57 events, suppressing the 
process (2) by a factor of 3 and making the 
expected background due to the process (3) 
very small, of the order of 10 events. Al- 
though, in contrast with the other processes, 
the Monte Carlo simulation of the process (3) 
is much less accurate due to nuclear inter- 
action of Kl, so the data analysis may not 
completely rely on this estimation. 

Characteristic feature of the process (1) 
is that each event must contain two pho- 
ton pairs with invariant masses of 77 and tt° 
mesons. Simulation shows that if the photons 
in an event are enumerated in descending or- 
der in energy, the photons from 77 — > 77 decay 
have the numbers of either 1 and 2 or 1 and 3. 
Corresponding experimental and simulated 
distributions in 777,12 and mi3 with an addi- 
tional requirement, that rest of photons con- 
tains a pair with \ra^ — m v o\ < 20 MeV/c 2 
and more stringent photon quality require- 
ment £ < —5. are shown in Fig.|l|. Back- 
ground estimations are based on PDG ta- 
ble value for <fi — > 777 decay branching ratio 
[ PHI , our measurements of e + e _ — > lutt cross 
section ]TT[| a nd cf> — > 7r°7r°7 decay branch- 
ing ratio ]l(|. Distribution of experimental 
events in 77742 + W13 shows an enhancement 
centered at 77-meson mass, while the sum of 
background processes is nearly flat in this re- 
gion. The sum of all simulated background 
processes, where each one was normalized to 
the number of events expected at a given inte- 
grated luminosity and total number </>-meson 
decays, describes the spectrum outside the 
enhancement quite well. If this excess of 
events is due to the decay (1), it corresponds 
to a 4> — > 777r°7 branching ratio of the order of 
7 ■ 10~ 5 , but it is difficult to extract accurate 
result from these inclusive spectra, because of 
poor signal to background ratio and unknown 
decay dynamics. 

Detailed study of the process (1) requires 
substantial suppression of background. It 
could be done using kinematic fit with in- 
termediate 77 and 7T° mesons taken into ac- 
count. For each event 777r°7 and 7T 7T 7 hy- 
potheses were tried and corresponding XjLo~j 
and x^-o T o 7 parameters calculated. To sup- 
press the processes (4) and (5) the follow- 
ing requirements were imposed: xf )7r o 7 < 20, 
X^-o^o., > 20, and for additional suppression 



the process (2): ( < —4. With these require- 
ments the contribution from the processes (4) 
and (5), which are relatively rare themselves, 
becomes negligible. 

The resulting spectra of the E jmax /E , en- 
ergy of the most energetic photon, are shown 
in Fig.^. Since the recoil photon in the pro- 
cess (2) has a narrow spectrum peaked at 360 
MeV, E imax /EQ must be more than 0.7 in 
this process. The spectrum in Fig.^Jb, where 
true 57 events are suppressed by an order 
of magnitude, is well reproduced by simula- 
tion of the process (2) alone, while the spec- 
trum in Fig.||a shows excess of experimen- 
tal events with C, < —4 over the simula- 
tion of (2). The background (2) still dom- 
inates at E imax / Eq > 0.7, so events with 
E~ ( max/E > 0.7 were rejected. From Fig.||a 
one can see, that this roughly halves the de- 
tection efficiency for the process (1). This 
also enhances dependence of the detection ef- 
ficiency on the invariant mass of rjir system, 
which varies from 1% at M vn o = 970 MeV/c 2 
up to 5% at M V7r o = 700 MeV/c 2 . 

The distributions of the events over pho- 
ton quality parameter £ are shown in FigJH 
Selection criteria here were the same as de- 
scribed above, except loosened requirement 
on photon quality: ( < 10. It can be 
seen, that simulation well describes the back- 
ground from the process (2), while the excess 
of events in Fig.||a at low £ is compatible 
with existence of the decay (1) with BR of 
the order of 10~ 4 . Also was checked the dis- 
tribution in y 2 • While the distributions 
of experimental events with < C < 10 is 
in a good agreement with simulation of the 
process (2), the distribution of the events 
with £ < —4 agrees with the simulation of 
the process (1) with a small contribution of 
the process (2). The total number of se- 
lected events satisfying the selection criteria 
is 25, from which 5 were estimated to be a 
background from the process (2). The corre- 
sponding branching ratio was calculated us- 
ing the r]Tr a invariant mass distribution of the 
selected events and detection efficiencies, ob- 
tained from simulation. The resulting value 
is (0.83 ± 0.23) • 10~ 4 , where error is a statis- 
tical one. For such stringent selection crite- 
ria the systematic error in detection efficiency 
could be as high as 15%, but at present level 
of experimental statistics the statistical error 
dominates. 

In order to study the rjir mass spectrum 
the E imax cut was dropped. This increased 
the background from the <fi — * VI decay, but 
made detection efficiency twice higher and 
less dependent on r/ir mass. The 777 back- 



ground was subtracted using events with M™ 
outside rj peak. The result is shown in Figjj. 
The invariant mass spectrum was fitted us- 
ing formulae from Ref.[Q and assuming that 
the decay is a pure — > a (980)7 transition. 
The result is also shown in Fig.|| It is seen, 
that although the errors are very large, the 
observed mass spectrum is consistent with 
007 decay mechanism. Constraint fit with 
9a n ir = 0.85g aK + K - H gives the following 
values of eto(980) meson parameters: 
M ao = 986 1 2 g MeV/c 2 , 
9 2 aK+K -/^=(1.5±0.5) GeV 2 . 

4 Conclusion 

In this work evidence of the new decay <p — » 
r]ir°j with a branching ratio of (0.83 ± 0.23) • 
10~ 4 is presented. The observed enhance- 
ment in the rjir invariant mass spectrum at 
large m v ^o is consistent with domination of 
ao7 decay mechanism, while relatively high 
decay rate agrees with the theoretical pre- 
dictions based on 4-quark model of ao(980)- 
meson §, §. 
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Figure 2: Spectra of the most energetic pho- 
ton in the event: a - photons quality ( < —4, 
b-0<£<10. Circles with error bars de- 
pict experimental data, shaded histogram - 
simulation of the process (2) and clear his- 
togram - simulated sum of (2) and (1) with 
a BR = 10" 4 . 



Figure 1 : Invariant masses of pairs of most 
energetic photons: circles with error bars - 
experimental data, a - simulated signal from 
4> — > ^7r°7 decay corresponding to a branch- 
ing ratio of 0.7 • 10~ 4 , b - estimated back- 
ground from the e + e~ — > urn and 4> ~ * 
VI ! fo (980)7 events, c - sum of a and b 




Figure 4: rjir invariant mass spectrum: 
points - experimental data, curve - optimal 
fit according to Ref. Q . 

Figure 3: Distribution of events over £: a - 
events with E imax / E$ < 0.7, b - events with 
0.7 < E imax /E < 0.8 circles with error bars 
- experimental data, shaded histogram de- 
picts simulation of the process (2) and clear 
histogram - simulated sum of (2) and (1) 
with &BR= 10~ 4 . 



